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We discuss a novel mechanism for the realization of supersolids in lattices in the absence of inter¬ 
site interactions that surprisingly works as well at unit filling. This mechanism, that we study for 
the case of the sawtooth lattice, is based on the existence of frustrated and un-frustrated plaquettes. 
For sufficiently large interactions and frustration the particles gather preferentially at un-frustrated 
plaquettes breaking spontaneously translational invariance, resulting in a supersolid. We show that 
for the sawtooth lattice the supersolid exists for a large region of parameters for densities above half- 
filling. Our results open a new feasible path for realizing supersolids in existing ultra-cold atomic 
gases in optical lattices without the need of long-range interactions. 

PACS numbers: 67.85.-d, 67.80.kb, 67.85.Hj 


I. INTRODUCTION 


Supersolids have attracted a large interest since they 
were proposecP^, due to their, apparently counter¬ 
intuitive, co-existence of both crystalline order and su¬ 
perfluidity. Superfluidity on top of the crystalline order is 
explained by the creation and delocalization of zero point 
defects (such as vacancies or interstitials) in a strongly 
interacting system. The search for the elusive supersolid 
remains a challenge 3 The claimed evidence of super¬ 
solidity in HeliunP, was subsequently explained by the 
shear modulus stiffening of solid 4 He213. 

Supersolidity may occur as well in lattices due to 
inter-site interactions. Lattice supersolidity has at¬ 
tracted an active theoretical interest as well 8 . Ultra-cold 
gases in optical lattices provide an interesting system for 
the realization of lattice supersolids. However, the re¬ 
quirement of sufficiently large inter-site interactions re¬ 
duces the possible scenarios for supersolidity to gases 
with long-range interactions, incl uding dipolar gased^ 
and Rydberg-dressed atomic gase*P- () 11 . These systems 
present however difficulties due to inelastic collisions in 
polar molecules 12 and short life times in Rydberg gases. 
An alternative to polar gases is provided by experiments 
with condensates in optical cavities, where infinitely long- 
range interactions between the condensed atoms are in¬ 
duced by two-photon processes. These interactions drive 
the Dicke phase transition that results in self-organized 
supersolids 13 . Although these experiments are realized in 
absence of optical lattices, the self-organized supersolids 
break spontaneously a discrete spatial symmetry, hence 
resembling the case of lattice super solids. 

Long-range interactions are however not necessary for 
the realization of lattice supersolids. Recent studies on 
frustrated lattices with flat bands, such as kagome lat¬ 
tices 14 and Creutz ladder^^, have discussed the possi¬ 
bility of observing supersolids without inter-site interac¬ 
tions. In the case of the Creutz ladder, supersolidity at 
incommensurate densities results from effective next-to- 
nearest neighbor hopping in the vicinity of the two fiat 
band region, and the doping of a Valence Bond Crystal. 

In this paper we discuss a novel mechanism that leads 



FIG. 1: (Color online) Sawtooth lattice discussed in the pa¬ 
per. The V-shaped plaquettes are un-frustrated, whereas the 
A-shaped plaquettes are frustrated, since we assume t > 0 
and t' < 0. The presence of frustrated and unfrustrated pla¬ 
quettes is crucial for the super solid mechanism discussed in 
the paper. 


to robust lattice supersolidity in the absence of inter-site 
interactions. The mechanism is based on the existence 
of frustrated and un-frustrated plaquettes. In this sense, 
although we study the specific case of a sawtooth lat¬ 
tice, we expect that supersolids may be also realized in 
other lattice geometries fulfilling this property. Super¬ 
solidity follows from the preferential occupation of un- 
frustrated plaquettes, which results in spontaneously bro¬ 
ken translational symmetry. Sawtooth lattices and other 
frustrated geometries may be realized using lattice shak- 
in P and similar laser arrangements as those recently 
employed for creating variable lattice geometries 18 20 . 
Hence our results open a new feasible path for realizing 
supersolidity in existing atomic gases without the need 
of long-range interactions. 


The structure of the paper is as follows. In Sec. |H| we 
introduce the sawtooth model studied in the paper. Sec¬ 


tion [Hi] discusses the non-interacting regime. Section |IV| 
is devoted to the ground state phase diagram, both at 
unit-filling and away from unit filling. In Sec. [V] we in¬ 
troduce a simple model that allows an intuitive under¬ 
standing of the super solid mechanism. Section [Vi] dis¬ 
cusses the roton instability responsible for the superffuid- 
to-supersolid transition. In Sec. |VII| we comment on the 
experimental realization of the sawtooth model, and on 
the signatures of the supersolid phase. Finally, we sum¬ 
marize our conclusions in Sec. Eml 
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FIG. 2: (Color online) Single particle dispersion for Model 0 
with U — 0 for different values of \t'\/t. 


II. MODEL 

We consider bosons in the sawtooth lattice as shown in 
Fig. B a), which is characterized by a hopping rate t' < 0 
along the lower leg, and a hopping rate t > 0 along the 
rungs. The change in sign of which may be achieved 
experimentally by e.g. lattice shakinjJ 17 ^ introduces ge¬ 
ometric frustration, a key ingredient below. The lower 
and upper legs have different coordination numbers, and 
hence constitute two distinct sublattices, which we de¬ 
note as A and B. For a sufficiently deep lattice, the 
system is described by the Bose-Hubbard model (BHM): 

n=-tY J (a\b i +b\a i+1 + H.c.) 

i 

+ t' y>K +1 + H.c.) + — y, n\(n\ — 1) (1) 

i v£{A,B},i 

where a\(b\) and a i (b/) are creation and annihilation op¬ 
erators for bosons at site i of leg A (B), and nf = a\a i 
('nf = b\b i ) is the number operator at site i of leg A (B). 
Contact-like interactions lead to the on-site interaction 
term, characterized by the coupling constant U. Hence¬ 
forth we consider all the physical quantities in units of 
t = 1 which makes them dimensionless. 


III. NON-INTERACTING REGIME 

At first, it is convenient to study the non-interacting 
regime (U = 0). In that case Model 0 may be di¬ 
agonalized in momentum (fc) space (note that due to 
the broken sublattice symmetry, k G [—7r/2 ,tt/2]): H = 

where the first (a) and sec¬ 
ond (/3) energy bands are characterized by the dispersion 
E a ,p(k) = — t' cos 2ky2t[cos 2 k + (t'/2t) 2 cos 2 2&] 1 / 2 , and 
the bosonic operators = cos Ok^k + sin Qkbk, and /3k = 
— sin OkCLk +cos Okbk, with tan 2 6k = (2t/1') cos k/ cos 2 k. 


The lowest band becomes flat at \t'\ = t/y /2 as shown 
in Fig. [2] When \t'\ < t/y/2 there is a single mini¬ 
mum in the lowest band, whereas two minima occur when 

\t'\>t/V2. 

As recently shown^, for the interacting case the flat 
band results in a solid phase at lattice filling p = 1/4, 
but the solid order breaks upon doping due to the prolif¬ 
eration of domain walls. In addition, at \t'\ = t/y /2 the 
minimum of E a (k ) changes from k = 0 to k = 7r/2. As 
a result, for non-interacting bosons, a transition occurs 
from a superfluid phase at k = 0 (SF 0 ) to a superfluid 
phase at k = 7r/2 (SF^)- Since a n / 2 = a n / 2 , in the 
SF ^/2 the B leg is depopulated. 

IV. GROUND-STATE PHASE DIAGRAM 

The ground state phase diagram is obtained by 
means of numerical density matrix renormalization 
group (DMRG) calculations 21 23 . For our calculation we 
consider up to L = 160 sites, with maximally 4 parti¬ 
cles per site, and retaining up to 500 states in the den¬ 
sity matrix. In order to obtain the signature of possible 
quantum phase transitions in the system we compute the 
corresponding order parameters. Particularly relevant is 
the density structure factor 

S(k) = ^y / e ik ^(n i n j ), ( 2 ) 

hi 

where ( riiUj ) is the density-density correlation between 
sites i and j. 

In order to quantify the critical point for the SFo to 
Mott-insulator (MI) transition discussed below, we calcu¬ 
late the Luttinger parameter A, which is the power-law 
exponent of the decay of the single particle correlation 
function T r = (a\aj). However, in order to obtain the 
value of K we evaluate the structure factor S(k) in the 
long wave-length limit {k -A 0), since in this limit 24 : 

S{k) oc K\k\/2'K. (3) 

We also compute the single particle gap to characterize 

the gapped MI phase as 

E G =y~Li~, ( 4 ) 

where p + (p~) is the chemical potential for adding (re¬ 
moving) one particle. Finally, the momentum distribu¬ 
tion function 

N{k) = ^y j e ik ^\a\a j ). (5) 

h3 

is particularly interesting as a possible experimental sig¬ 
nature of the supersolid phase discussed below. 

In the following subsections we discuss the ground- 
state phase diagram at unit filling, and away from unit 
filling. 
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FIG. 3: (Color online) Phase diagram of Model 0 for p — 1 
as a function of \t'\/t and U/t (see text). 


A. Unit filling 

Figure [3] shows the ground state phase diagram of 
model 0 at unit filling, p = 1. As discussed above, 
in the non-interacting limit the system exhibits two gap¬ 
less superfluid phases, SFo and SF^/ 2 - In addition to 
these phases, the system opens as expected a MI phase 
at small C7, which at the flat-band point is maintained 
all the way down to vanishing U . The second additional 
phase is, in contrast, rather unexpected. As we discussed 
below, this phase is a supersolid (SS), which rather sur¬ 
prisingly occurs at p = 1 and in the absence of inter-site 
interactions. 

In the following we discuss in detail the different phases 
for the case U/t = 3. In Fig. [4] we plot the thermody¬ 
namic limit values of the single-particle gap 0 • The MI 
region is characterized by a finite gap, whereas Eg —> 0 



FIG. 4: (Color online) Extrapolated gap Eg as a function of 
\t'\/t. The inset shows the exponential to power-law decay 
of the single particle correlation function F r at the MI-SS 
boundary. Only odd r is shown to avoid oscillations. 


(a) (b) 




FIG. 5: (Color online) (a) Extrapolated values of K (black 
squares) with respect to \t'\/t for U/t = 3. The K — 2 
line (red dashed) is drawn to find the critical point for the SFo- 
MI transition; (b) Scaled momentum distribution N(k)L~ 3 U 
for different lengths as a function of \t'\/t for U/t — 3. 


in all other regions. The smooth gap opening at small 
\t'\/t marks the Berezinskii-Kosterlitz-Thouless (BKT) 
SFo-to-MI transition, characterized by a Luttinger pa¬ 
rameter K = By performing a finite size scaling of 
K computed using Eq. (|3| , we obtain that the transition 
for SF 0 to MI phase occurs at \t'\/t ~ 0.1 (see Fig.[5ja)). 
The BKT transition is further confirmed by the scaling 
of the momentum distribution, N(k = 0) oc L l ~ 27 . 
In Fig. [HJb) we show that the values of N(k = 0)L -3 / 4 
for different lengths (L=80, 120, 160) intersect at the 
transition point \t'\/t ~ 0.1. 

In contrast, the sudden vanishing of Eq for larger 
\t'\/t (at \t'\/t = 1.4 in Fig.|4| does not match with a BKT 
transition. The gapped to gapless transition is, how¬ 
ever, confirmed by the single-particle correlation function 
T r = (ajai+ r ), which decays exponentially in the MI and 
algebraically for \t'\/t > 1.4 (inset of Fig. |4|. This tran¬ 
sition is likely to be a weak first order transition which 
could not be confirmed in our numerics. 

The superfluid region opening at the large \t'\/t side 



FIG. 6: (Color online) Number fluctuation between odd and 
even sites of the A and B legs, the inset shows A n and A d- 
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FIG. 7: (Color online) (Top) Scaled structure factor S(k = 
7t/2) with respect to \t'\/t for U/t = 3.0; (Bottom) Finite-size 
scaling of S(k = 7t/2) for U/t = 3.0 and different \t'\/t values. 



FIG. 8: (Color online) Phase diagram of the sawtooth chain 
as a function of p and \t'\/t for U/t — 2.0. Black lines with 
squares show the gapped regions. The SS phase is bounded by 
the red circles. The SF ^/2 phase occupies a large portion of 
the phase diagram below the SS region and for \t'\/t > l/\/2, 
also around p — 1. The region bounded by magenta triangles 
shows a macroscopic jump in the density. The brown dashed 
lines are two representative cuts at \t'\/t = l/y/2 and 1.5 (see 
text). The red dot corresponds to the flat-band point where 

|F|/* = 1/V2 


k = ±7t/2, confirmed by extrapolation to the thermody¬ 
namic limit. In Fig. [7] (top) we show, for U/t = 3, the 
extrapolated value of S(k = n/ 2), which is finite in the 
region 1.4 < \t'\/t < 3.1, i.e. the SS phase. The extrapo¬ 
lation of 5(fc = 7r/2) is shown in detail in Fig. [7] (bottom). 


of the MI turns out to be a SS phase. Figure [6] shows 
the average density avg^ dd even (avg^ ddeven ) at even and 
odd sites in the A (B) leg. The difference of densities 
at both legs is not surprising due to the asymmetry be¬ 
tween the legs. More interesting is the behavior within 
a given leg. Whereas for avg^ dd = avg^ ven for all \t'\/t, 

aV 9odd + av 9even in the SS re S ion ( L4 < \?\/t < 3.1 
in Fig. [6]) indicating a clear density modulation in this 
region. The spontaneously broken translational symme¬ 
try along the A leg is characterized by the amplitude 
A n = avg^ dd — avg^ ven of the odd-even modulation (inset 
of Fig. [6|. The broken symmetry translates as well in a 
difference of bond-kinetic energy Ad = avg^ dd —avg^ ven , 

where avg^ dd even = \ E ie A odd , eve M-i ai + a l bi + 1 + 

JT.c.), which is finite in the SS region (inset of Fig. [6|. 
The finite A d relates with the V-type dimerization dis¬ 
cussed below. The end of the SS region and the onset of 
the SF ^/2 phase is marked by the vanishing of both A n 
and A£>, and the depopulation of the B leg. 

The SS phase can be further confirmed by a finite peak 
in S(k) at non-zero wave vector k. However, we note 
that the sawtooth lattice breaks a trivial translational 
symmetry due to which S(k) shows a peak at k = ±7r. 
Therefore, the presence of the SS phase in the system 
is revealed by a finite peak in the structure factor at 


B. Away from unit filling 

At this point we consider incommensurate filling. This 
discussion is particularly relevant, since the presence of 
an overall harmonic trap, typical of experiments on atoms 
in optical lattices, results in an inhomogeneous density 
distribution. In Fig. [8] we depict our results for U/t = 2. 
In addition to the gapped density-wave phase at p = 0.25 
at the flat-band pointwe obtain gapped phases at 
p = n/4, for all integer n > 1 (black lines and squares 
in Fig. [8|. The gapped phases at these densities are ex¬ 
pected to occur at the large interaction limit. The phases 
at pm (2n + l)/4 with n > 1 are similar to the phase at 
p = 1/4, being comparatively narrow due to small inter¬ 
action. On the contrary, the gapped phase at p = 1/2 is a 
trivial insulator due to a filled band. Even more relevant 
is the robustness of the SS phase, which extends over a 
large parameter region, for densities p > 0.5. 

We illustrate the behavior of the system for the case 
of two representative values, \t'\ = l/y/2 and 1.5 (brown 
dashed lines in Fig. |8|. 

In Fig. [9j we plot p as a function of the chemical po¬ 
tential /i for \t'\/t = l/y/2. Due to band flatness at that 
point, at a critical g there exists a macroscopic jump in 
the density up to p = 0.25. This macroscopic jump exists 
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FIG. 9: (Color online) p(p) for L = 40, U/t — 2, and \t'\/t = 
l/y/2 (blue curve) and 1.5 (red curve). 


for a small region around \t'\/t = l/\/2. This region is 
marked in the phase diagram by magenta triangles. At 
p = 0.25 a gap appears marked by a plateau in p(/x). 
Then, the system becomes a gapless and compressible 
superfluid. Further increase in p induces additional gaps 
as discussed above. In between the gapped plateaus the 
system is a superfluid. 

The dependence of p(p) is very different for \t'\/t = 1.5, 
since for that value the band is not flat. As a result the 
system starts at p 0 in the SF ^/2 phase. In Fig. [9] 
we observe that the curve p(p) is continuous with two 
kinks at p ~ 0.8 and 1.2. These kinks correspond to the 
transition to the SS phase. The abrupt growth of p in the 
SS phase, shows that the SS phase is highly compressible. 

Finally, Fig. 10 shows S(k = 7r/2) for \t'\/t = l/y/2 
and 1.5. For |rj ~[t = 1.5, S(k = 7r/2) increases in the 
SS phase, 0.8 < p < 1.2. In contrast, for \t'\/t = 1/v^, 
S(k = 7r/2) remains small except at three peaks at p = 
0.25,0.75 and 1.25. The two small peaks at p — 0.75 
and 1.25 are due to the presence of two small density- 



FIG. 10: (Color online) Structure factor S(k = 7t/2) for dif¬ 
ferent fillings p for L = 40, U/t = 2, and \t'\/t » l/y/2 (blue 
circles) and 1.5 (red squares). 


wave phases shown as black squares in Fig. [8] The sharp 
peak at p = 0.25 corresponds to the density-wave phase 
reported in Reff® 


V. WEAK COUPLING LIMIT 

An intuitive insight on the emergence of the SS is ob¬ 
tained from the classical limit of Model 0 . We assume 
each site to be in a coherent state with a well defined 
density and phase, 0 or 7r, corresponding to the two pos¬ 
sible minima of E a (k). We consider a simplified model 
in which bj = ip a 2 j = £, and a> 2 j+i = X, such that we 
allow for both a possible density imbalance between the 
A and B legs, and for an even-odd asymmetry in the A 
leg. We may hence minimize the energy, which without 
loss of generality may be calculated for a four-site unit 
cell: 


(H) = -4 t(£x + XV) ~ 4 

+ U{? + 2x 4 + r/ 4 )-m(£ 2 + 2x 2 +?? 2 ) (6) 

Within this approach the phase diagram splits into three 
regions (Fig. fll^d)). For \t'\/t < \/\[2 the three 
coefficients have the same sign, and the particles oc¬ 
cupy both A and B sites corresponding to the SFq- 


phase (Fig. 11'a)). For small U/t and \t'\/t > 1/V% the 


B sites depopulate (p = 0), the density is homogeneous 
in the A sites (|x| = |£|), and sigri(x) ^ sign(£), corre¬ 
sponding to the SF ^/2 phase (Fig. flTlb)). A sufficiently 
strong repulsive interaction U/t > ( U/t) c redistributes 
population to the B sites. However, how the particles re¬ 
distribute in the A leg is crucially determined by the ex¬ 
istence of frustrated and un-frustrated plaquettes in the 
sawtooth lattice. In order to minimize kinetic energy par¬ 
ticles favor the un-frustrated V-shaped plaquettes of the 
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FIG. 11: (Color online) Sketch of the SFo (a), SFtt /2 (b) 
and SS (c). The circle size is proportional to the local den¬ 
sity, whereas d= denote the sign, (d) Phase diagram in the 
weak-coupling limit as function Up and \t'\/t for a fixed (but 
arbitrary) filling. The SF 0 to the SS phase is shown as a 
solid line. The dotted and dashed lines depict, respectively, 
the SFtt/ 2-SS transition line obtained from the variational 
approach § and the roton instability. At t — 0 SS-order 
vanishes (thick solid line) resulting in an highly degenerate 
ground-state, since particles may occupy the uncoupled B 
sites in an arbitrary configuration. 
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FIG. 12: (Color online)Densities of the classical model for 
pU — 2 and (a) t = 1 and t' — —0.5 and (b) t — 0.8 and 
t! — — 1 as a function of the chemical potential p. The dotted 
(blue) line depicts the average density of the B sites, the solid 
(black) and dashed (red) lines show densities of odd and even 
A sites. 


sawtooth forming V-shaped dimers. As a result particles 
break the translational symmetry spontaneously, prefer¬ 
ably occupying every second V-plaquette (Fig. [TT] (c)), 
which leads to a density modulation in the A sites, that 
characterizes, as mentioned above, the SS phase. Note 
that this simple picture also predicts a finite Ap in the 
SS region, as observed in the numerics. 

As discussed above, already the classical or weak - 
coupling limit gives an intuitive insight into the formation 
of the SS phase. We provide in this section some addi¬ 
tional details. In Fig. [12] we show typical densities as 
function of /i, obtain using Model |6|. While in the SFo 
and SF ^/2 phases the A-sites do not exhibit any density 
modulation, when entering the SS-phase the avg^ ven and 
av 9^dd curves separate from each other. The amplitude 
of the modulation on the A sites in the limit p £, \t'\, U 
is given by £ — x = is suppressed with increasing 

repulsion U. Within this approach the population of the 
B sites in the SF ^/2 phase is strictly zero. 

We determine the SF 0 to SS transition as first-order, 
whereas the SF ^/2 to SS transition is a second-order 
phase transition exhibiting a discontinuity in in second 
derivative of the ground-state energy. 


VI. ROTON INSTABILITY. 

Insight on the nature of the SF^/ 2 -SS transition is ob¬ 
tained in the limit in which U <C Ep(k) — E a {k ) for all 
k. In that case we may project Model M onto the lowest 
energy band: 

H~Y,E a (k)ala k + ^'EtiP’ k '- 9a 

k q,k,k' 

(7) 

with flit = nticos(^) + Note that, 

although the on-site interactions are contact-like, the ef¬ 
fective interactions are momentum dependent. 



FIG. 13: (Color online)Bogoliubov spectrum e(k) for \t'\/t = 
1 and (from top to bottom) pU — 0, 0.1, 0.2, 0.3, 0.4, 0.5, 
0.5278. At k — 0 a roton-like minimum develops, which 
touches zero at pU ~ 0.5278. 


Starting from Eq. 0 and after expanding the in¬ 
teraction part up to second order terms assuming con¬ 
densation at q = 7t/ 2, i.e. a q ~ v+ a q the total 
Hamiltonian may be written up to constants as H ~ 

T,k>o A ( k ) («4A + & -k & -k) +B(k) (a{a f _ k + a k a- k ) 
with 


A(k) = E a (k ) + pU [cos20 q (cos20k — cos 2 6 q ) 

+ (1 + COS 2 0 q COS 26k)] 

B(k) = pU (1 + cos 26 q cos 26k) (8) 


This Hamiltonian may be readily diagonalized us¬ 
ing a Bogoliubov transformation (3k = cosh — 
sinh 7 fcQ:i fe yielding the Bogoliubov spectrum of excita¬ 
tions, e(k) 2 = (^E^ik) + 2Up(l + cos20*;)^ E A {k) with 
E^(k) = E A (k) - 2 1' + U p(cos 26 k - 1). 

In Fig.[l3]we depict the corresponding energy spectrum 
emerging for finite pU. As usual the spectrum exhibits 
a linear (phonon-like) dispersion for k close to n/2. In¬ 
terestingly, for finite pU it acquires a local minimum at 
k = 0, that resembles the roton dispersion minimum of 
superfluid He^, and that occurs, as for dipolar conden¬ 
sate^, due to the momentum dependence of the inter¬ 
actions. For a critical value of pU the roton-like min¬ 
imum reaches zero energy, becoming unstable, marking 
the transition to the SS. As shown in Fig. [TT] (d) the criti¬ 
cal pU for roton instability agrees well with the SF^/ 2 -SS 
transition line obtained by the classical model. Hence, we 
can conclude that the SF ^/2 is destabilized through roton 
instability that leads to the SS phase. 
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FIG. 16: (Color online) Sawtooth lattice formed by an in¬ 
coherent superposition between a Kagome lattice 20 and an 
additional lattice V(f) — sin 2 (^^-ky^j + \ sin 2 (^Q-ky^ with 
k the laser wavenumber. Darker regions mean lower potential. 


FIG. 14: (Color online)Momentum distribution N(k) in dif¬ 
ferent regions of the phase diagram shown in Fig.[3]for U — 3. 


VII. EXPERIMENTAL SIGNATURE 


In the section we briefly discuss about the possible sig¬ 
natures of the SS phase in optical lattice experients us¬ 
ing ultracold atoms. The sawtooth lattice can be created 
by suitably using a superlattice potential on top of the 
Kagome lattice which has been created recentl}^ (see 
Fig. 16). 


In order to obtain the signature of the SS phase we 
compute the momentum distribution N(k) using Eq. |5j 
As expected, in SFo (SF^/ 2 ) N(k) has a peak at k = 
0 (k = 7T/2). However, an intriguing feature appears in 
the SS phase which shows, interestingly, peaks at k = 0 
and k = tt/2 as shown in Fig. [l4j We also plot the peak 
strengths as a function of \t'\/tlit k = 0 and k = 7t/2 in 
Fig. [15] which shows the transition from SFq - MI - SS - 



FIG. 15: (Color online) iV(0) (red squares) and N( 7t/2) (blue 
circles) in different regions of the phase diagram shown in 
Fig. [3] at 17 = 3. 


SF ^/2 for U = 3 and p = 1. Hence, the appearance of 
the SS phase may be directly monitored in time-of-flight 
experiments from the multi-peaked momentum distribu¬ 
tion. 


VIII. CONCLUSIONS 


We have discussed a novel mechanism for the forma¬ 
tion of lattice supersolids for the particular case of a saw¬ 
tooth lattice. The mechanism is based on the selective 
population of un-frustrated plaquettes in the presence of 
frustrated and unfrustrated plaquettes, and hence we ex¬ 
pect supersolids in other lattices fulfilling that property. 
We have shown that the supersolid exists for a broad 
range of lattice fillings, including various commensurate 
fillings, in particular unit filling. Interestingly the super¬ 
solid may be revealed not only by in-situ measurements, 
but by monitoring the momentum distribution in time- 
of-flight measurements. Since frustrated lattices, in par¬ 
ticular sawtooth, may be realized using state of the art 
techniques, our results open hence a new feasible path 
for realizing supersolids in existing experiments with ul¬ 
tra cold atoms in optical lattices, without the need of 
long-range interactions. 
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